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Abstract 



_ We present a new approach to enhancing Answer Set Programming (ASP) with Con- 

fi I straint Processing techniques which allows for solving interesting Constraint Satisfaction 

Problems in ASP. We show how constraints on finite domains can be decomposed into 
logic programs such that unit-propagation achieves arc, bound or range consistency. Ex- 
periments with our encodings demonstrate their computational impact. 
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^D ' 1 Introduction 

O' 

Answer Set Programming (ASP; IBaral 2003P has been put forward as a power- 
ful paradigm to solve Constraint Satisfaction Problems (CSP) in (jNiemela 1999^ . 
Indeed, ASP has been shown to be a useful in various applications, among them 
planning (|Lifschitz 1999^ . model checking ( [Heljanko and Niemela 2003P , and bio- 
C^ ' informatics (jBaral et al. 2004|) . and decision support for NASA shuttle controllers 

( [Nogueira et al. 2001[ ). It combines an expressive but simple modelling language 
with high-performance solving capacities. In fact, modern ASP solvers, such as 
ciasp (jGebser et al. 2007a|) . compete with the best Boolean Satisfiability (SAT; 
IBiere et al. 2009|) solvers. An empirical comparison of the performance of ASP and 
traditional Constraint Logic Programming (CLP: [JafFar and Maher 1994^ on solv- 



ing CSP conducted by Dovier et al. (2005) shows ASP encodings to be more com 



pact, more declarative, and highly competitive, but also revealed shortcomings: 
non-Boolean constructs, like resources or functions over finite domains, in par- 
ticular global constraints, are more naturally modelled and efficiently handled by 
Constraint Processing (CP: IDechter 20031 IRossi et al. 2006[) systems. 

This led to the integration of CP techniques into ASP. Similar to Satisfiability 
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Modulo Theories fSMT: lNieuwenhuis et al. 2006p . the key idea of an integrative ap- 
proach is to incorporate theory-specific predicates into propositional formulas, and 
extending an ASP solver's decision engine for a more high-level proof procedure. 
Recent work on combining ASP with CP was conducted in (jBaselice et al. 20051 
IMeUark od et al. 20081 iMellarkod and Gelfond 200"8t and ()Gebser et al. 2009p . While 
[McUark od and Gelfondl both view ASP and CP solvers as blackboxes. lGebser et al.l 
embed a CP solver into an ASP solver adding support for advanced backjumping 



and conflict-driven learning techniques. Balduccini (2009) and Jarvisalo et al. (2009) 



cut ties to ad-hoc ASP and CP solvers, and principally support global constraints. 



Dal Palu et al. (2009) put further emphasis on handling constraint variables with 
large domains, and presented a strategy which only consider parts of the model that 
actively contribute in supporting constraint answer sets. However, each system has 
a subset of the following limitations: either they are tied to particular ASP and CP 
solvers, or the support for global constraints is limited, or communication between 
the ASP and CP solver is restricted. 

This paper introduces a translational approach to Constraint Answer Set Solving 
rather than an integrative one. Motivated by the success of SAT-based constraint 
solvers, such as the award- winning system Sugar (jTamura et al. 2006^ . we show 
how to enhance ASP with Constraint Processing techniques through translation to 



ASP. A first study was conducted in Gebser et al. (2009) with the system xpanda 



for representing multi-valued propositions in ASP. One of the key contributions 
of our work is an investigation of constraint decomposition techniques in the new 
field of Constraint Answer Set Programming, illustrated on the popular all-different 
constraint. The resulting approach has been implemented in the new preprocessor 
inca. Empirical evaluation demonstrates its computational potential. 

The remainder of this paper is organized as follows. We start by giving the back- 
ground notions of ASP and Constraint Satisfaction. Various generic ASP encodings 
of constraints on finite domains and proofs of their properties are given in Section |31 
In Section|4l we empirically evaluate our approach and compare to existing research. 
Section [5] draws conclusions. 

2 Background 

2.1 Answer Set Programming 

A (normal) logic program over a set of primitive propositions ^ is a finite set of 
rules of the form 

h i~ ai,. . .,am, not «,„+!, . . . , not a„ 

where < m < n and h, ai (£ A is an atom for 1 < i < n. A literal a is an 
atom a or its default negation not a. For a rule r, let head(r) — h he the head of 
r and body{r) = {oi, . . . , «„, not a^n+i, ■ ■ ■ , not a„} the body of r. Furthermore, 
define hody{r)^ = {ai, . . . , «,„} and body(r)^ — {a„i+i, ■ ■ ■ , an}. The set of atoms 
occurring in a logic program 11 is denoted by atom(n), and the set of bodies in 11 
is body{Il) = {body{r) \ r G H}. For regrouping bodies sharing the same head a, 
define body{a) = {body{r) | r e 11, header) = a}. 
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The semantics of a program is given by its answer sets. A set X C yl is an 
answer set of a logic program 11 over A, if X is the C-minimal model of the reduct 
(|Gelfond and Lifschitz 19881) 

n-^ = {head{r) ■h- body{r)+ | r £ H, body{ry n X = 0}- 

The semantics of important extensions to normal logic programs, such as choice 
rules, integrity and cardinality constraints, is given through program transforma- 
tions that introduce additional propositions (cf. ISimons et al. 2002]) . A choice rule 
allows for the non-deterministic choice over atoms in {hi, . . . ,hk} and has the fol- 
lowing form: 

{/ii, . . . , /ifc} <- ai, . . . , a™, not a,n+i,- . . ,not a„ 

An integrity constraint 

<- Oi, . . . , flm, not flm+i, ...,not a„ 

is an abbreviation for a rule with an unsatisfiable head, and thus forbids its body 
to be satisfied in any answer set. A cardinality constraint 

■It- k{di, . . . , a„} 

is interpreted as no answer set satisfies k literals of the set {ai, . . . , a„}. It can 
be transformed into (Jj integrity constraints r such that body{r) C {ai, . . . , a„} 
and \body{r)\ = k. ISimons et al.l provide a transformation that needs just 0(nk) 
rules, introducing atoms l{ai,j) to represent the fact that at least j of the literals 
with index > j, i.e. the literals in {di, . . . , a„}, are in a particular answer set can- 
didate. Then, the cardinality constraint can be encoded by an integrity constraint 
■<— l{di,k) and the three following rules, where 1 < i < n and 1 < j < k: 

l{di,j) ^ l{d^+ij) /(a,,j + 1) ^ a,,/(a,+i,j) l{d^,l) ^ di 

Notice that both transformations are modular. Alternatively, modern ASP solvers 
also incorporate propagators for cardinality constraints that run in 0{n). 



2.2 Nogoods 

We want to view inferences in ASP as unit-propagation on nogoods. Following 



Gebser et al. (2007b), inferences in ASP rely on atoms and program rules, which 
can be expressed by using atoms and bodies. Thus, for a program 11, the domain 
of Boolean assignments A is fixed to dom{A) = atom (11) U body(n). 

Formally, a Boolean assignment A is a set {cri, . . . , cr„} of signed literals ai 
for 1 < i < n of the form Ta or Fa where a G dom{A). Ta expresses that 
a is assigned true and Fa that it is false in A. (We omit the attribute Boolean 
for assignments whenever clear from the context.) The complement of a signed 
literal a is denoted by a, that is Ta = Fa and Fa = Ta. In the context of ASP, 
a nogood (jDechter 2003P is a set J = {ci, . . . , cr„} of signed literals, expressing a 
constraint violated by any assignment A such that (5 C A. For a nogood i5, a signed 
literal a & 5, and an assignment A, we say that S is unit and a is unit-resulting if 
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Input: A set V of nogoods, and an assignment A. 

Output: An extended assignment, and a status (either conflict or success) 



repeat 

if (5 C A for some 5 £ V then 

return (A, conflict); 
T.^{SeV\S\A = {a},a ^ A}; 
if E / then let a G <5 \ A for some 5 e E in 
A <- AU (a); 
until E = 0; 
return (A, success); 



Fig. 1. The unit-propagation algorithm. 

S\ A — {a}. Let A'^ = {a G dom{A) \ Ta e ^4} the set of true propositions and 
A^ = {a e dom{A) \ Fa £ A} the set of false propositions. A total assignment, 
that is A'^ U A^ — dom{A) and A^ U A'^ = 0, is a solution for a set A of nogoods 
if (5 g A for all S e A. 



As shown in Lee (20051, the answer sets of a logic program H correspond to 
the models of the completion of H that satisfy the loop formulas of all non-empty 
subsets of atom{lV). For j3 — {ai, . . . , a„j, not Um+i, ■ ■ ■ , not a„} G body{Il), define 

_ f {Tai,. . . ,Ta,„,Fa,„+i,. . .Fa„,F/3}, 
^ ~ \ {Fai, T/3}, . . . , {Fa™, T/3}, {Ta„+i, T/3}, . . . , {Ta„, T/3} 

Intuitively, the nogoods in A^ enforce the truth of body /? iff all its literals are 
satisfied. For an atom a G atom{Il) with body{a) = {/3i, . . . , /?fc}, let 

r {F/3i,...,F/3fc,Ta}, 
"^ \ {T/3i,Fa},...,{T/3fc,Fa} 

Then, the solutions for An = U/3e body (n) ^/^ U {Jaeatorn{n)^a correspond to the 
models of the completion of H. Loop formulas, expressed in the set of nogoods An, 
have to be added to establish full correspondence to the answer sets of H. Typically, 
solutions for An U An are computed by applying ConHict-Driven Nogood Learning 
(CDNL: lGebser et al. 2007b[) . This combines search and propagation by recursively 
assigning the value of a proposition and using unit-propagation (Fig.[T]) to determine 
logical consequences of an assignment (jMitchell 2005p . 

Example 1 

Consider the set of nogoods V — {{Tai, Fa2, Taa, Ta4}, {Fai, Ta4}, {Fa3,Ta4}} 

and the assignment A = {Ta4}. Unit-propagation extends A by {Tai, Ta2, Taa}. 



2.3 Constraint Satisfaction and Consistency 

A Constraint Satisfaction Problem is a triple ( V, D, C) where F is a set of variables 
V = {vi, . . . , Vn}, D is a, set of finite domains D — {Di, . . . , _D„} such that each 
variable Vi has an associated domain dom{vi) = Di, and C is a set of constraints. 



Following Rossi et al. (2006), a constraint c is a pair {Rs,S) where Rs is a fc-ary 
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relation on the variables in 5 C V^ , called the scope of c. In other words, Rs is 
a subset of the Cartesian product of the domains of the variables in S. To access 
the relation and the scope of c define range{c) — Rs and scope{c) — S. For a 
(constraint variable) assignment A : V —^ [Jvev dom,{v) and a constraint c — 
{Rs, S) with S = {vi, . . . , Ufc), define A{S) — {A{vi), . . . , A{vk)), and call c satisfied 
if A{S) G range{c). A binary constraint c has \scope{c)\ — 2. For example, vi ^ V2 
ensures that vi and V2 take different values. A global (or n-ary) constraint c has 
parametrized scope. For example, the all-different constraint ensures that a set of 
variables, {vi, . . . , w„} take all different values. This can be decomposed into 0{n^) 
binary constraints, Vi ^ Vj for i < j. However, as we shall see, such decomposition 
can hinder inference. An assignment ^ is a solution for a CSP iff it satisfies all 
constraints in C. 

Constraint solvers typically use backtracking search to explore the space of par- 
tial assignments, and prune it by applying propagation algorithms that enforce a 
local consistency property on the constraints after each assignment. A binary con- 
straint c is called arc consistent iff when a variable vi G scope{c) is assigned any 
value di e dom{vi), there exists a consistent value d2 G dom{v2) for the other vari- 
able V2. An n-ary constraint c is hyper-arc consistent or domain consistent iff when 
a variable Vi G scope{c) is assigned any value di G dom{vi), there exist compatible 
values in the domains of all the other variables dj G dom(vj) for all i<j<n, j^i 
such that (di, . . . , (i„) G range{c). 

The concepts of bound and range consistency are defined for constraints on or- 
dered intervals. Let min{Di) and m,ax{Di) be the minimum value and maximum 
value of the domain Di . A constraint c is bound consistent iff when a variable Vi is 
assigned di G {min(dom(vi)) , max{dom{vi))} (i.e. the minimum or maximum value 
in its domain), there exist compatible values between the minimum and maximum 
domain value for all the other variables in the scope of the constraint. Such an 
assignment is called a bound support. A constraint is range consistent iff when a 
variable is assigned any value in its domain, there exists a bound support. Notice 
that range consistency is in between domain and bound consistency, where domain 
consistency is the strongest of the four formalisms. 



2.4 Constraint Answer Set Programming 



Following Gebser et al. (2009), a constraint logic program 11 is defined as logic 
programs over an extended alphabet distinguishing regular and constraint atoms, 
denoted by A and C, respectively, such that head{r) G A for each r G H. Constraint 
atoms are identified with constraints via a function 7 : C — > C, and furthermore, 
define 7(6") = {7(c) | c G C"} for C C C. For a (constraint variable) assign- 
ment A define the set of constraints satisfied by A as satc{A) = {c \ A{scope{c)) G 
range{c), c G C}, and the constraint reduct as 

n^ = {head{r) <- body{r)\j( | r G 11, 

-l{hody{r)+\c) C satc{A), -i{hody{r)-\c) n satc{A) ^%}- 
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Then, a set X C ^ is a constraint answer set of 11 with respect to ^, if X is an 
answer set of 11"^ . 

In the translational approach to Constraint Answer Set Solving, a constraint logic 
program is compiled into a (normal) logic program by adding an ASP decompo- 
sition of all constraints comprised in the constraint logic program. The constraint 
answer sets can then be obtained by applying the same algorithms as for calculat- 
ing answer sets, e.g. CDNL. Since all variables will be shared between constraints, 
nogood learning techniques as in CDNL exploit constraint interdependencies. This 
can improve propagation between constraints. 

3 Encoding Constraint Answer Set Programs 

In this section we explain how to translate constraint logic programs with multi- 
valued propositions into a (normal) logic program. There are a number of choices of 
how to encode constraints on multi- valued propositions, e.g. a constraint variable v, 
taking values out of a pre-defined finite domain, dom{v). In what follows, we assume 
dom{v) = [1, d] for all u G V to save the reader from multiple superscripts. 

3.1 Direct Encoding 

A popular choice is called the direct encoding (jWalsh 2000p . In the direct encoding, 
a propositional variable e(v, i), representing v = i, is introduced for each value i 
that can be assigned to the constraint variable v. Intuitively, the proposition e{v, i) 
is true if v takes the value i, and false if v takes a value different from i. For 
each V, the truth- assignments of atoms e{v,i) are encoded by a choice rule (1). 
Furthermore, there is an integrity constraint (2) to ensure that v takes at least one 
value, and a cardinality constraint (3) that ensures that v takes at most one value. 

(1) {e(w,l),...,e(«,rf)}^ 

(2) <— not e{v, 1), . . . , not e{v, d) 

(3) ^2{e{v,l),...,e{v,d)} 

In the direct encoding, each forbidden combination of values in a constraint is 
expressed by an integrity constraint. On the other hand, when a relation is rep- 
resented by allowed combinations of values, all forbidden combinations have to be 
deduced and translated to integrity constraints. Unfortunately, the direct encoding 
of constraints hinders propagation: 

Theorem 1 (| Walsh 2000\ 

Enforcing arc consistency on the binary decomposition of the original constraint 

prunes more values from the variables domain than unit-propagation on its direct 

encoding. 

3.2 Support Encoding 

The support encoding has been proposed to tackle this weakness (jGent 2002| . A 
support for a constraint variable v to take the value i across a constraint c is the 
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set of values {ii, . . . , «„} C dom(v') of another variable in v' £ scope{c) \{v} which 
allow V — i, and can be encoded as follows, extending (1-3): 

<— e{v^ i), not e{v' , ii), . . . , not e{v\ i,n) 

This integrity constraint can be read as whenever v — i, then at least one of its 
supports must hold. In the support encoding, for each constraint c there is one 
support for each pair of distinct variables v, v' € scope(c), and for each value i. 

Theorem 2 \Gent 2002} 

Unit-propagation on the support encoding enforces arc consistency on the binary 

decomposition of the original constraint. 

We illustrate this approach on an encoding of the global all-different constraint. For 
variables w, v' and value i it is defined by the following 0{n^d) integrity constraints: 

<— e(?;, i), not e{v\ 1), . . . , not e{v' , i — 1), not e{v' , i + 1), . . . , not e{v' , d) 
To keep the encoding small, we make use of the following equivalence (e) 

e{v , i) = not e{v , 1), . . . , not e{v , i ^ 1), not e{v , j + 1), . . . , not e{v , d) 
covered by (2-3) and get 

■^ e{v, i), e{v' , i)- 

Observe, that this is also the direct encoding of the binary decomposition of the 
global all-different constraint. However, this observation does not hold in general 
for all constraints. As discussed in the Background section of this paper, we can 
express above condition as 0{d) cardinality constraints: 

(4) i-2{e{vi,i),...,e{vn,i)} 

Corollary 1 

Unit-propagation on (1-4) enforces arc consistency on the binary decomposition of 

the global all-different constraint in 0{nd^) down any branch of the search tree. 

Proof 

From the definition of cardinality constraints, (4) ensure that for all distinct v,v'e 
scope{c) any value i is not taken by both v and v' . These integrity constraints 
correspond to the support encoding of the global all-different constraint since (2-3) 
cover the equivalence (e). By Theorem[2l unit-propagation on this support encoding 
enforces arc consistency on the binary decomposition of the all-different constraint. 
For each of the n variables, there are 0{d) nogoods resulting from (1-3) that 
can be woken 0{d) times down any branch of the search tree. Each propagation 
requires 0(1) time. Rules (1-3) therefore take 0{nd^) down any branch of the 
search to propagate. There are 0{nd) nogoods resulting from (4) that each take 
0(1) time to propagate down any branch of the search tree. The total running time 
is given by 0{nd'^) -H 0{nd) = 0{nd'^). D 
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3.3 Range Encoding 

In the range encodings a prepositional variable r(w, I, u) is introduced for all [I, u] C 
[1, d\ to represent whether the value of v is between I and u. For each range [Z, u], 
the following 0{nd?) rules encode v G [Z, m] whenever it is safe to assume that 

V ^ [1, / — 1] and V ^ [m + 1, d], and enforce a consistent set of ranges such that 

V e[l,u]^ V <^[l -l,u]Av e[l,u + 1]: 

(5) r(i), Z, u) ^— noi r{v, 1, Z — 1), not r{v, m + 1, d) 

(6) -^ r(w, Z — 1, u), not r{v, /, m) 

(7) •<— r(w, Z, u + 1), not r{v, /, m) 

Constraints are encoded into integrity constraints representing conflict regions. 
When the combination vi G [/i, mi], . . . , w„ G [?«,««] violates the constraint, the 
following rule is added: 

<- r(ui,/i,Mi),...,r(w„,/„,u„) 

Theorem 3 

Unit-propagation on the range encoding enforces range consistency on the original 

constraint. 



Proof 

Suppose we have a set of ranges on the domains of the constraint variables in which 
no unit-propagation is possible and no domain is empty. Consider any constraint 
variables ui , . . . , u^ , . . . , w„ and value di such that there is no bound support in 
Vi,. . . , Vi-i, Wi+i, ■ ■ . Vn for Vi = di, i.e. there are no compatible values for the other 
variables Vj distinct from Vi where Vj G [Ij, Uj]. Hence, all instantiations such that 
Vi — di are in a conflict region vi G [h, ui] C [l[, u[], . . . , Vn € [In, w„] C [/^, u^]. For 
each Vj we have Tr{vj,lj,Uj) representing Vj G [lj,Uj\. Then, the binary nogoods 
{Tr(wj, Ij, Uj),Fr{vj, /, — 1, Uj)} and {Tr(w,-, Ij, Uj),Fr{vj, Ij, Uj + l)} resulting from 
(6) and (7) are unit, and eventually we get Tr{vj, l!j, u'A for [Ij, Uj] C [/', u'A. But 
then the nogood {Tr{vi, l[, u[), . . . ,Tr(w„, Z^, u'n)} encoding the conflict region is 
unit and forces ¥r{vi, l[, u[) representing Vi ^ [l[, u[]. By nogoods resulting from (6) 
and (7) we get di is not in the domain of Vi , and the domains are bound consistent 
as required. Since at least one value must be in each domain, encoded in (5), we 
have a set of non-empty domains which are range consistent. D 

A propagator for the global all-different constraint that enforces range consistency 



pruning Hall intervals has been proposed in Leconte (1996 1 and encoded to SAT in 
[Bessiere et al. (2009| ). An interval [/, u] is a Hall interval iff Iji" | dom{v) C [I, u]}\ = 
u — I + 1. In other words, a Hall interval of size k completely contains the domains 
of k variables. Observe that in any bound support, the variables whose domains are 
contained in the Hall interval consume all values within the Hall interval, whilst 
any other variable must find their support outside the Hall interval. 



A Translational Approach to Constraint Answer Set Solving 9 

Example 2 

Consider the global all- different constraint over the variables {wi, ^27 ^3, ^^4} with 
dom{vi) = {2,3}, dom{v2) = {1,2,4}, dom{vs,) = {2,3}, dom{v4) = {1,2,3,4}. 
[2, 3] is a Hall interval of size 2 as the domain of 2 variables, vi and U3, is completely 
contained in it. Therefore we can remove [2, 3] from the domains of all the other 
variables. This leaves V2 and v^ with a domain containing values 1 and 4. 

The following decomposition of the global all-different constraint will permit us to 
achieve range consistency via unit propagation. It ensures that no interval [I, u] can 
contain more variables than its size. 

(8) ^ u- 1 + 2 {r{vi, ;, m), . . . , r(w„, I, u)} 

This simple decomposition can simulate a complex propagation algorithm like 
ILecontel s with a similar overall complexity of reasoning. 

Corollary 2 

Unit-propagation on (5-8) enforces range consistency on the global all-different 

constraint in 0{nd^) down any branch of the search tree. 

Proof 

Clearly, the cardinality constraints (8) reflect all conflict regions such that no Hall 
interval [Z, u] can contain u — I -\-2 variables, that are more variables than its size. 
Hence, (8) is a range encoding of the global all-different constraint. By Theorem |31 
unit-propagation on this encoding enforces range consistency on the global all- 
different constraint. 

There are 0{nd^) nogoods resulting from (5-7) that can be woken 0{d) times 
down any branch of the search tree. Each propagation requires 0(1) time. Rules 
(5-7) therefore take 0{nd^) down any branch of the search to propagate. There are 
0{nd^) nogoods resulting from (8) that each take 0{l) time to propagate down any 
branch of the search tree. The total running time is given by 0{nd^) + 0{nd^) = 

o{nd^). a 



3.4 Bound Encoding 

A last encoding is called the bound encoding (jCrawford and Baker 1994]) . In the 
bound encoding, a propositional variable b{v,i) is introduced for each value i to 
represent that the value of v is bounded by i. That is, v < i if T6(w, i), and v > i 
if Fb{v,i). Similar to the direct encoding, for each v, the truth-assignments of 
atoms b{v, i) are encoded by a choice rule (9). In order to ensure that assignments 
represent a consistent set of bounds, the condition v<i^v<i + lis posted as 
integrity constraints (10). Another integrity constraint (11) encodes v < d, that at 
least one value must be assigned to v: 

(9) {b{v,l),...,b{v,d)}i- 

(10) ^ b{v,i),not b{v,i + l) Vi £ [1, d - 1] 

(11) -^notbiv^d) 
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Constraints are encoded into integrity constraints representing conflict regions sim- 
ilar to the range encoding. When aU combinations in the region 

h < Vi < Ui, . . . , In < Vn < Un 

violate a constraint, the following rule is added: 

^ b{vi, Ml), . . . , b{vn, Un), uot b{vi, k), . . . , not 6(w„, In) 

Theorem 4 

Unit-propagation on the bound encoding enforces bound consistency on the original 

constraint. 



Proof 

Suppose we have a set of bounds on the domains of the constraint variables in which 
no unit-propagation is possible and no domain is empty. Consider any constraint 
variable Vi such that if Vi is assigned its minimum domain value l^ + l or its maximum 
domain value Ui there are no compatible values of the other constraint variables 
vi, . . . , Wi_i, Wi+i, ■ ■ ■ ,Vn between their minimum h+l, . . . , 4_i+l, Zj+i+l, . . . , 4i+l 
and their maximum domain values ui, . . . , Wi^i, Mi+i, . . . ,Un, respectively. First, we 
analyse the case Vi = Ui, that is, all instantiations such that Vi = Ui are in a conflict 
region I'l < h < vi < ui < u'l, . . . , l'^ < In < Vn '^ Un < u'^. For each Vj we have 
F 6 ( w, , Ij ) and T 6 ( w, , Mj ) , representing Vj > Ij and Vj < Uj . Then, the binary nogoods 
{Tb{vj,lj — l),Fb{vj,lj)} and {Tb{vj,Uj),Fb{vj,Uj + 1)} resulting from (10) are 
unit, and eventually we get Fb{vj, I') for L' < Ij as well as T6(w,, u') for m' > Uj. 
But then the nogood {F6(wi, ^{), T&(wi, m(), . . . ,F6(?;„, Z4),T&(w„, u^)} encoding 
the conflict region is unit and forces Tb{vi, l[) representing v^ < l[. Since l[ < Ui 
and by the nogoods resulting from (10) we get Ui is not in the domain of Vi. 

The second case, where Vi is assigned its minimum domain value k + 1, is sym- 
metric, and we conclude that the domains are bound consistent as required. Since 
at least one value must be in each domain, resulting from (11), we have a set of 
non-empty domains which are bound consistent. D 

In order to get a representation of the global all-different constraint that can only 
prune bounds, the bound encoding for variables is linked to (8) as follows: 

(12) r{v, I, u) <— not b{v, I — 1), b{v, u) 

(13) ^ r{v,l,u),b{v,l-l) 

(14) <— r{v,l,u),not b{v,u) 

Corollary 3 

Unit-propagation on (8-14) enforces bound consistency on the global all-different 

constraint in 0{n<P) down any branch of the search tree. 
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Proof 

This result follows from Corollary [5] and Theorem S] Observe that the decomposi- 
tions for range and bound consistency both encode the same conflict regions. 

For each of the n variables, there are 0{d) nogoods resulting from (9-11) that 
can be woken 0{d) times down any branch of the search tree. Each propagation 
requires 0(1) time. Rules (9-11) therefore take 0{n(P) down any branch of the 
search to propagate. Furthermore, there are 0{nd?) nogoods resulting from (8) 
and (12-14) that each take 0(1) time to propagate down any branch of the search 
tree. The total running time is given by 0{nd?). D 

Note that an upper bound h can be posted on the size of Hall intervals. The resulting 
encoding with only those cardinality constraints (5) for which u— l + l < h detects 
Hall intervals of size at most /i, and therefore enforces a weaker level of consistency. 

4 Experiments 

To evaluate our decompositions, we conducted experiments on encodingqll of CSP 
containing all-different and permutation constraints. The global permutation con- 
straint is a special case of all-different when the number of variables is equal to the 
number of all their possible values. A decomposition of permutation extends (4) by 

<— not e{vi, «),..., not e(u„, i) 

or (8) by the following rule where 1 < I < u < k: 

■^ d — u + I {not r{vi, I, u), . . . , not r(u„, /, u)} 

This can increase propagation. Our translational approach to Constraint Answer Set 
Solving has been implemented within the prototypical preprocessor inca. Although 
our semantics is propositional, inca compiles constraint logic programs with first- 
order variables, function symbols, and aggregates, etc. in linear time and space, 
such that the logic program can be obtained by a grounding process. Experiments 
consider inca^ in different settings using different decompositions. We denote the 
support encoding of the global constraints by S, the bound encoding of the global 
constraints by i?, and the range encoding of the global constraints by R. To explore 
the impact of small Hall intervals, we also tried Bk and Rk , an encoding of the global 
constraints with only those cardinality constraints (8) for which u ^ I + 1 < k. The 
consistency achieved by Bk and Rk is therefore weaker than full bound and range 
consistency, respectively. 

We also include the integrated systems clingcon^ (0.1.2), and ezcsjo (1.6.9) in 
our empirical analysis. While clingcon extends the ASP system clingo^ (2.0.2) with 
the generic constraint solver gecodqj (2.2.0), ezcsp combines the grounder gringo^ 
(2.0.3) and ASP solver clasp^ (1.3.0) with sicstuci (4.0.8) as a constraint solver. 

^ http://potassco.sourceforge.net/ 

^ http : //krlab . cs . ttu . edu/'marcy/ezcsp/ 

^ http://www.gecode.org/ 

* http://www.sics.se/sicstus/ 
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Table 1. Data on time and space for selected translations. 









S 






B 




R 




n 


time 


atoms 


rules 


atoms 


rules 


atoms 


rules 




12 


0.1 


236 


151 


1,829 


1,907 


1,061 


2,423 


PHP 


14 


0.1 


304 


177 


2,876 


2,981 


1,630 


3,877 




16 


0.1 


380 


203 


4,263 


4,399 


2,375 


5,823 




8 


0.3 


655 


36,525 


1,297 


41,565 


3,217 


38,237 


QG5 


10 


0.9 


1,219 


109,066 


2,421 


118,946 


7,121 


111,446 




12 


2.3 


2,039 


267,643 


4,057 


284,755 


13,849 


270,067 




4 


0.6 


564 


74,945 


941 


78,848 


4,135 


81,174 


DWn 


6 


3.0 


1,130 


363,115 


1,983 


373,002 


12,581 


381,724 




8 


9.2 


1,888 


1,122,549 


3,409 


1,142,132 


28,355 


1,163,810 



Note that dingo stands for clasp on gringo and combines both systems in a mono- 
lithic way. Since inca is a pure preprocessor, we select the ASP system dingo (2.0.3) 
as its backend to provide a representative comparison with dingcon and ezcsp. The 
behaviour of xpanda is simulated by setting S and, therefore, is not considered in 
our study. We also do not separate time spend on grounding and solving the prob- 
lem, since the grounder's share of the overall runtime is generally insignificant on 
our benchmarks. 

To compare the performance of Constraint Answer Set solvers against pure CP 
systems, we also report results of gecode (3.2.0). Its heuristic for variable selection 
was set to a smallest domain as in dingcon. All experiments were run on a 2.00 GHz 
PC under Linux. We report results in seconds, where each run was limited to 600 
s time and 1 GB RAM. 

Space Complexity. Data on the size of selected translations shown in Table [T] con- 
firms our theoretical results. For PHP and DWn instances (description follows), 
the number of atoms in the support (bound/range) encoding is bounded by 0{n^) 
(0(n3)). The number of rules is ©(n) (©(n^)) for PHP, 0{n^) (©(n^)) for DWn 
due to constraints represented in the direct encoding. An n x n table is modelled in 
QG5, raising the number of atoms to 0{n^) {0{n'^)), and rules to 0{n*) {0{n*)). 



4-1 Pigeon Hole Problems 

The Pigeon Hole Problem (PHP) is to show that it is impossible to put n pigeons 
into n — 1 holes if each pigeon must be put into a distinct hole. Clearly, our bound 
and range decompositions are faster compared to weaker encodings (see Table [2]). 
On such problems, detecting large Hall intervals is essential. 



4-2 Quasigroup Completion 

A quasigroup is an algebraic structure (Q,-)j where Q is a set and ■ is a binary 
operation on Q such that for every pair of elements a,b £ Q there exist unique 



A Translational Approach to Constraint Answer Set Solving 13 

Table 2. Runtime results in seconds for PHP. 

n S Bi B2 B3 B R3 R ezcsp clingcon gecode 



10 


5.4 


0.7 


0.5 


0.1 


0.0 


0.2 


0.0 


1.8 


1.4 


0.9 


11 


46.5 


3.5 


1.5 


1.0 


0.0 


1.9 


0.0 


16.7 


15.2 


9.0 


12 


105.0 


14.8 


7.1 


3.9 


0.0 


2.6 


0.1 


183.9 


172.5 


104.1 


13 


— 


91.4 


68.6 


25.4 


0.1 


30.4 


0.0 


— 


— 


— 


14 


— 


— 


350.1 


125.0 


0.0 


196.9 


0.1 


— 


— 


— 


15 


— 


— 


— 


— 


0.1 


— 


0.1 


— 


— 


— 


16 


— 


— 


— 


— 


0.1 


— 


0.1 


— 


— 


— 



Table 3. Average times over 100 runs on QCP. Timeouts are given in parenthesis. 
% S Bs B Rs R ezcsp clingcon gecode gecodesc 



10 


2.6 


5.0 


8.2 


6.0 


7.3 


29.6 (7) 


9.7 (4) 


2.2 (4) 


0.5 (1) 


20 


2.4 


5.0 


8.0 


6.2 


7.2 


21.3 (20) 


6.2 (5) 


5.0 (4) 


0.9 (3) 


30 


2.3 


4.8 


7.9 


6.1 


7.1 


10.3 (30) 


12.9 (13) 


2.9 (13) 


1.1 (5) 


35 


2.3 


4.8 


7.9 


6.1 


7.0 


21.6 (24) 


11.2 (17) 


14.1 (13) 


6.2 (7) 


40 


2.3 


4.7 


7.8 


6.0 


6.9 


51.6 (29) 


23.1 (22) 


11.7 (20) 


5.7 (9) 


45 


2.3 


4.7 


7.8 


5.9 


6.8 


36.3 (35) 


14.7 (28) 


17.7 (25) 


6.3 (13) 


50 


2.3 


4.6 


7.7 


5.9 


6.8 


36.1 (50) 


21.2 (37) 


25.1 (32) 


6.3 (18) 


55 


2.3 


4.5 


7.6 


5.8 


6.7 


61.4 (51) 


24.4 (44) 


19.6 (41) 


30.9 (29) 


60 


2.2 


4.4 


7.5 


5.6 


6.6 


60.2 (63) 


31.4 (56) 


36.0 (51) 


27.2 (35) 


70 


2.2 


4.2 


7.1 


5.1 


6.0 


70.0 (66) 


30.2 (50) 


28.0 (45) 


17.0 (27) 


80 


2.1 


4.0 


6.7 


4.7 


5.5 


16.2 (18) 


4.2 (18) 


17.2 (13) 


7.0 (7) 


90 


2.1 


4.0 


6.7 


4.7 


5.5 


1.4 


2.6 (1) 


0.4 (1) 


3.2 



elements x,y ^ Q which solve the equations a ■ x = b and y ■ a — b. The or- 
der n of a quasigroup is defined by the number of elements in Q. A quasigroup 
can be represented by an n x n-multiplication table, where for each pair a, b the 
table gives the result of a • b, and it defines a Latin square. This means that each 
element of Q occurs exactly once in each row and each column of the table. The 
Quasigroup Completion Problem (QCP) is to determine whether a partially filled 
table can be completed in such a way that a multiplication table of a quasigroup is 
obtained. Randomly generated QCP has been proposed as a benchmark domain for 



CP systems by Gomes and Selman (19971 since it combines the features of purely 



random problems and highly structured problems. Table |3] compares the runtime 
for solving QCP problems of size 20 x 20 where the first column gives the percent- 
age of preassigned values. We included gecode with algorithms that enforce bound 
and domain consistency, denoted as gecode bc a-nd gecode oc (not shown), in the 
experiments. Our analysis exhibits phase transition behaviour of the systems ezcsp., 
clingcon, gecode, and gecode bc, while our Boolean encodings and gecode dc solve 
all problems within seconds. Interestingly, learning constraint interdependencies as 
in our approach is sufficient to tackle QCP. In fact, most of the time for S, Bk, Rk 
is spent on grounding, but not for solving the actual problem. 
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Table 4. Runtime results in seconds for QEP. 





n 


S 


Bi 


B:i 


B 


i?3 


R 


ezcsp 


clingcoD 


gecode 




7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.6 


65.0 


189.8 


0.6 


QGl 


8 


19.0 


5.9 


4.7 


19.8 


6.4 


4.7 


— 


— 


— 




9 


— 


139.4 


152.0 


234.6 


27.6 


466.9 


— 


— 


— 




7 


1.7 


1.7 


1.7 


1.8 


1.7 


1.8 


46.1 


1.5 


1.2 


QG2 


8 


46.6 


9.6 


10.6 


37.7 


11.7 


14.8 


— 


— 


— 




9 


— 


246.0 


55.7 


88.3 


119.7 


213.4 


— 


— 


— 




7 


0.2 


0.2 


0.2 


0.3 


0.2 


0.3 


3.2 


1.0 


0.0 


QG3 


8 


0.4 


0.4 


0.5 


0.5 


0.5 


0.5 


4.3 


9.0 


0.2 




9 


10.2 


7.4 


9.5 


16.5 


11.0 


12.8 


— 


— 


18.2 




7 


0.2 


0.2 


0.2 


0.3 


0.3 


0.3 


2.8 


0.7 


0.1 


QG4 


8 


0.5 


0.6 


0.7 


0.9 


0.8 


0.7 


27.9 


36.8 


0.3 




9 


1.3 


1.0 


2.1 


3.0 


1.1 


0.9 


442.1 


288.8 


3.7 




8 


0.4 


0.4 


0.4 


0.5 


0.4 


0.4 


6.9 


5.3 


0.0 




9 


0.7 


0.8 


0.8 


0.9 


0.8 


0.8 


249.2 


— 


0.0 


QG5 


10 


1.6 


1.5 


1.6 


1.9 


1.6 


1.6 


— 


— 


0.2 




11 


2.1 


2.2 


2.4 


3.4 


2.8 


2.4 


— 


— 


0.8 




12 


27.0 


6.2 


9.1 


12.4 


8.4 


10.4 


— 


— 


16.4 




8 


0.4 


0.4 


0.5 


0.5 


0.5 


0.4 


0.8 


— 


0.0 




9 


0.7 


0.7 


0.8 


0.9 


0.8 


0.8 


1.2 


— 


0.0 


QG6 


10 


1.2 


1.4 


1.5 


1.8 


1.6 


1.5 


10.5 


— 


0.1 




11 


2.7 


2.8 


4.0 


4.2 


3.9 


4.8 


125.5 


— 


1.2 




12 


32.0 


12.9 


25.6 


36.4 


25.7 


50.6 


— 


— 


24.6 




8 


0.4 


0.4 


0.4 


0.6 


0.5 


0.5 


1.1 


— 


0.1 


QG7 


9 


0.7 


1.0 


1.2 


1.7 


1.2 


1.4 


9.1 


— 


0.9 




10 


6.7 


3.2 


5.2 


8.0 


4.7 


4.6 


— 


— 


22.0 



4-3 Quasigroup Existence 

The Quasigroup Existence Problem (QEP) is to determine the existence of cer- 



tain interesting classes of quasigroups. We follow Fujita et al. (1993) and look at 
problems QGl to QG7 that were target to open questions in finite mathematics. 
We represent them in the direct encoding which weakens the overall consistency. 
Furthermore, we add the axiom a ■ n > a — 1 where n is the order of the desired 
quasigroup, to avoid some symmetries in search space. We also assume quasigroups 
to be idempotent, that means a ■ a — a. QEP has been proposed as a benchmark 



domain for CP systems in Gent and Walsh (1999). All axioms have been modelled 



in ezcsp and gecode using constructive disjunction, and in S, B^, Rk and clingcon 
using integrity constraints. Table |4] demonstrates that both constructive disjunc- 
tion and integrity constraints have a similar behaviour, as for ezcsp and clingcon 
on benchmark classes QGl to QG4. On harder instances, conflict-driven learning 
appears to be too costly for clingcon. Additional experiments revealed that clingcon 
without learning performs like ezcsp. On the other hand, our decompositions benefit 
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Table 5. Runtime results in seconds for GGP. 



DWu 


S 


Bi 


53 


B 


Rz 


R 


ezcsp 


clingcon 


gecode 


3 


11.4 


3.8 


5.7 


8.7 


6.0 


10.4 


6.5 


66.9 


1.8 


4 


1.3 


2.0 


1.5 


3.2 


3.0 


2.5 


0.6 


0.1 


0.1 


5 


4.5 


5.0 


4.5 


13.5 


12.5 


31.4 


1.0 


2.0 


0.1 


6 


7.2 


11.0 


17.6 


47.7 


21.3 


110.2 


1.2 


— 


7.2 


7 


23.8 


28.3 


67.9 


227.9 


60.0 


432.9 


18.0 


— 


— 


8 


48.4 


68.4 


— 


207.8 


58.4 


356.8 


4.3 


— 


— 


9 


82.8 


106.5 


200.4 


486.6 


227.4 


— 


390.5 


— 


— 



from learning constraint interdependencies, resulting in runtimes that outperform 
all other systems including gecode on the hardest problems. 

4-4 Graceful Graphs 

A labelling / of the nodes of a graph {V,E) is graceful if / assigns a unique 
label f{v) from {0,1, . . . ,\E\} to each node v £ V such that, when each edge 
{v, w) £ E is assigned the label \f{v) — /(w)|, the resulting edge labels are distinct. 
The problem of determining the existence of a graceful labelling of a graph (GGP) 



has been modelled as a CSP in Petrie and Smith (20031. Our experiments consider 



double- wheel graphs DWn composed by two copies of a cycle with n vertices, each 
connected to a central hub. Table [5] shows that our encodings compete with ezcsp 
and outperform the other comparable systems, where the support encoding per- 
forms better than bound and range encodings. In most cases, the branching heuris- 
tic used in our approach appears to be misled by the extra variables introduced in 
Bk and Rk- That explains some of the variability in the runtimes. 

5 Conclusions 

We have provided a new translation-based approach to incorporating Constraint 
Processing into Answer Set Programming. In particular, we investigated various 
generic ASP decompositions for constraints on finite domains and proved which 
level of consistency unit-propagation achieves on them. Our techniques were for- 
mulated as preprocessing and can be applied to any ASP system without chang- 
ing its source code, which allows for programmers to select the solvers that best 
fit their needs. We have empirically evaluated their performance on benchmarks 
from CP and found them outperforming integrated Constraint Answer Set Pro- 
gramming systems as well as pure CP solvers. As a key advantage of our novel 
approach we identified CDNL, exploiting constraint interdependencies which can 
improve propagation between constraints. Future work concerns a comparison to 
INiemelal s encoding (jNiemela 1999| lYou and Hou 2004)) . and encodings of further 
global constraints useful in Constraint Answer Set Programming. 

Acknowledgements. We are grateful to Martin Gebser and Torsten Schaub for useful 
discussions on the subject of this paper. 
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